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We have calculated the properties of SrTiOs (STO) using a wide array of density functionals 
ranging from standard semi-local functionals to modern range-separated hybrids, combined with 
several basis sets of varying size/quality. We show how these combination's predictive ability varies 
significantly, both for STO's cubic and antiferrodistortive (AFD) phases, with the greatest variation 
in functional/basis set efficacy seen in modeling the AFD phase. The screened hybrid functionals we 
utilized predict the structural properties of both phases in very good agreement with experiment, 
especially if used with large (but still computationally tractable) basis sets. The most accurate 
results presented in this study, namely those from HSE06/modified-def2-TZVP, stand as the most 
accurate modeling of STO to date when compared to the literature; these results agree well with 
experimental structural and electronic properties as well as providing insight into the band structure 
alteration during the phase transition. 
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I. INTRODUCTION 

Strontium titanate (SrTiOs; STO) is a complex ox- 
ide perovskite of great technological interest for its 
superconductivity,^ blue-light emission,^ photovoltaic 
effect,^ and so on. Under normal conditions, bulk 
SrTiOs crystallizes in a cubic perovskite structure; it sub- 
sequently undergoes a second order phase transition at 
Tc=105 K to a tetragonal structure with slightly rotated 
oxygens around the z-axis, known as the antiferrodis- 
tortive (AFD) phase (see Fig. 1). Many of the interest- 
ing properties of STO, either in bulk or in superlattices 
formed with other metal oxides, are believed to be caused 
by the cubic to AFD phase transition. Examples of this 
attribution are STO's superlattice's high Tc supercon- 
ductivity^^ and its colossal magnetoresistivity.^ First- 
principles calculations (see Ref. 8 and references therein) 
have indicated that the strain-induced competition be- 
tween octahedral rotation modes and the lattice distor- 
tion in metal oxide superlattices are behind these inter- 
esting properties. Thus, there is a considerable need^'^^ 
for precise theoretical calculations of the structural and 
electronic properties of complex oxides, as well as accu- 
rate estimation of the phase transition order parameters, 
to understand and eventually exploit these phenomena. 

The phase transition of STO is governed by two order 
parameters. The primary order parameter is the rotation 
angle of the TiOe octahedra {0). The experimentally 
measured^^ octahedral rotation of AFD STO is 1.4° at 
77 K and increases as the temperature drops toward the 
maximum measured value of 2.1° at 4.2 K. The octahe- 
dron's rotation is believed to be almost complete^ ^ at 
around 50 K, where ^=2.01±0.07° was reported. The 
secondary order parameter is the tetragonality of the unit 
cell (c/a), which increases from 1.00056^^ to 1.0009^^ as 
the temperature decreases from 65 K to 10 K.^^ The AFD 
phase can also appear in thin films of STO^^^^ at much 




FIG. 1. (Color online) SrTiOs unit cells for the a) cubic 
phase and b) antiferrodistortive phase; b) shows the TiOe 
octahedra's rotation around the [001] axis. Sr atom are in 
green, Ti are in blue and O are in red. The Oi (equatorial) 
and O2 (axial) labels denote the non-rotating and rotating 
oxygens, respectively. 



higher Tc than the bulk, depending on the substrate 
used, the thickness of deposited STO film, the strain and 
the lattice mismatch. For example, 10 nm of STO de- 
posited on LaAlOs (LAO) undergoes a transition to the 
AFD phase at ^ 332 K. 

As the simplest metal oxide perovskite, STO has been 
extensively studied in the last decades with different ah 
initio schemes. ^^^^ However, it is still a challenging ma- 
terial for theory; only a few of the previously published 
works have been able to accurately describe the struc- 
tural and electronic properties of the both phases of STO. 
The balance of this section will consist of a brief review 
of the theoretical work performed to date. 

Sai and Vanderbilt^^ carried out one of the first LDA 
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calculations on STO using a plane-wave basis and ultra- 
soft pseudopotentials. LDA predicted an exaggerated 
tetragonal AFD phase of STO, with octahedral rotation 
angles of 6°, significantly overestimating the 2.1° rotation 
measured experimentally.^^ Using LDA with other basis 
sets^^ shows similar issues, predicting rotations up 8.4°. 

Wahl et al.^^ used a plane- wave basis while simulat- 
ing STO with LDA26, PBE27'28 and PBEsoL^^^^o ^g^^ 
Section II for further descriptions of these density func- 
tionals). LDA underestimated experimental lattice con- 
stants, while PBE overestimated them; both methods 
had band-gaps that were seriously underestimated com- 
pared to experiment. This underestimation is well known 
for these functionals; see e.g. Ref. 31 and references 
therein. PBEsol was found to reproduce accurately the 
experimental structure, but considerably underestimated 
the band gaps. For the AFD phase, the octahedral an- 
gle was found to be very sensitive to the functional 
used; all three overestimate the AFD deformation, with 
LDA worse than PBE and PBEsol splitting the differ- 
ence. Rondinelli et al.^^ applied the LSDA-ht/ correc- 
tion to cubic STO and found that while it corrects the 
band gap, the calculated octahedral rotation angle re- 
mains overestimated at 5.7°. To date, none of the post- 
DFT corrections which benefit band-gaps have success- 
fully corrected the octahedral rotation overestimation, 
and many authors attribute this to the argument pro- 
posed by Sai and Vanderbilt^^ stating that this can be 
caused by the exchange and correlation terms in DFT 
not capturing quantum point fluctuations. 

Piskunov et al.^^ conducted one of the most complete 
and comprehensive ah initio studies of STO, using Gaus- 
sian basis sets specifically optimized for modeling STO 
crystals. This study of STO showed problems when 
modeling with pure DFT or pure HF, namely underesti- 
mated and overestimated band gaps, respectively; this is 
a well known problem. Hybrid functionals, specifically 
B3PW^^ and BSLYP,^^ gave more reasonable results, 
with direct band gaps overestimated by 5% for B3PW 
and 3.5% for B3LYP compared to experiment and indi- 
rect band gaps overestimated by 12% for B3PW and 10% 
for B3LYP. (We will demonstrate that an important part 
of this overestimation can be attributed to the basis set 
employed; see section III.) The hybrid functionals also 
gave the best agreement with experiment for the lattice 
constant and the bulk modulus, and generally did better 
than semilocal functionals in all categories. This success 
of hybrid functionals motivated more detailed calcula- 
tions^^"^^ of the properties of the cubic and AFD phases 
of STO, again using the optimized basis set of Piskunov 
et al.^^ and the B3PW functional. 

Next, Wahl et al.^^ applied the Heyd-Scuseria- 
Ernzerhof^'''^^ screened Coulomb hybrid density func- 
tional (HSE) in a plane- wave basis set. HSE performed 
exceptionally well, doing much better than any of semilo- 
cal functionals, as it gave a very accurate estimate of 
both the structural and electronic properties of the cubic 
phase. HSE also showed excellent agreement with the 



experimental octahedral angle and tetragonality of the 
unit cell which constitute, to our knowledge, the most 
accurately computed STO properties available in the lit- 
erature for both phases, prior to the current study. 

As noted above, hybrid functionals have proved their 
effectiveness in studying metal oxides, but they are com- 
putationally much more demanding than semilocal func- 
tionals. While it would be ideal to do high accuracy ab 
initio on metal oxide superlattices using complete basis 
sets and large supercells, this is prohibitively expensive 
at the current level of computer power. Screened hy- 
brid functionals with only short range exact exchange 
are computationally less demanding; they allow the use 
of large supercells, especially when used with localized 
basis sets such as Gaussian functions. We hope to use 
the most effective methods/basis sets from this study on 
more complicated metal oxide systems, and thus we have 
concentrated on methods and basis sets that would be 
practical for those systems as well as the systems cur- 
rently under consideration. 

This paper focuses on two tightly linked problems. We 
are interested in the degree of completeness (or size) of 
the localized basis set necessary to correctly simulate 
both phases of STO, and in the efficacy of recently de- 
veloped functionals (including screened hybrids) in pre- 
dicting the properties of STO. To discuss these issues, 
the paper proceeds as follows: In Section II, we briefly 
describe the technical details before turning in Section 
III to the basis set optimization/modification technique 
we used to make standard basis sets PBC-compatible. In 
Section IV, we report the results of semilocal and range 
separated hybrid functionals applied to the cubic and the 
AFD phases of STO. We show also how the quality of ba- 
sis set affected the accurate prediction of the octahedral 
rotation angle in the AFD phase of STO. Finally, we dis- 
cuss the results of our best functional/basis set combina- 
tion for STO, comparing them with previously published 
theoretical and experimental data, with special emphasis 
on the effect of varying the range separation parameter 
in the screened functionals. 



II. COMPUTATIONAL DETAILS 

All calculations were performed using a develop- 
ment version of the Gaussian suite of programs, 
with the periodic boundary condition (PBC)^^ code 
used throughout. A wide array of functionals were 
applied, including: the Local Spin Density Approxi- 
mation^^ (LSDA), the generalized gradient approxima- 
tion (GGA) corrected functional of Perdew, Burke and 
Ernzerhof^^'^^ (PBE), the reparametrization of PBE 
for solids, PBEsol,^^'^^ the revised meta-GGA of Tao, 
Perdew, Staroverov and Scuseria^^'^^ (revTPSS), and fi- 
nally a modern and highly parametrized meta-GGA func- 
tional, M06L.^^'^^ Two screened hybrid functionals were 
also tested, namely the short-range exact exchange func- 
tional of Heyd, Scuseria, and Ernzerhof^^'^^ (HSE, with 
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the 2006 errata, also referred to as HSE06) and the ex- 
act exchange in middle-range functional^^ of Henderson, 
Izmaylov, Scuseria, and Savin (HISS).^^'^^ Because reg- 
ular hybrids with unscreened exact exchange like B3LYP 
and B3PW have higher computational cost compared to 
screened hybrids we decided to exclude them from this 
test. 

Gaussian basis sets of different quality have been tested 
for their ability to simulate the properties of STO; the 
details of these tests and the modification of the basis set 
was detailed enough to merit its own section, section III. 

A few numerical considerations should be mentioned 
here. During the initial (or exploratory) calculations for 
the AFD phase, we found some dependence of octahe- 
dral rotation angle {0) on initial atomic positions. After 
further investigation, this can be attributed to the geo- 
metric optimization convergence criteria. Since is so 
small, very stringent convergence criteria is required. 
Another modified (versus the default) setting was that 
a pruned integration grid for DFT of (99, 590) was em- 
ployed, which corresponds to the Gaussian option "ultra- 
fine". Note that this grid is big enough for this system 
to avoid any of the instabilities with M06L reported in 
the literature with small grids. ^^'^^ To ensure this, we 
tested M06L with a larger grid, without noticing any 
modification in the calculated properties. Thus while "ul- 
trafine" is sometimes insufl&cient for M06L, it is not for 
this system. Other numerical settings in GAUSSIAN were 
left at the default values, e.g. integral cut-offs, k-point 
meshes^^, SCF convergence criterion, and the like. 

Finally, the geometry of each phase is worth discussing 
briefiy. The starting configuration for the cubic phase 
(see Figure 1(a)) consisted of the perovskite primitive 
cell containing 5 atoms at the experimental lattice con- 
stants^ (ao= 3.890 A). For the AFD phase, we couldn't 
simply use the 5 atom tetragonal unit cell with rotated 
oxygens and the lattice parameters set to a = 6 7^ c. 
A 20 atoms supercell/simulation cell was necessary, as 
the phase transition requires a rotation of every pair of 
neighboring TiOe octahedra in opposite directions (fig- 
ure 1(b)). Thus, the volume of the AFD supercell is 
about four times the volume of the cubic phase with 
tetragonal lattice constants a* = 6* = >/2a and c* = 2c, 
with a and c being the lattice parameters of the 5 atoms 
tetragonal unit cell in the AFD phase. The starting AFD 
structure of STO was taken from the experimental struc- 
ture of Jauch et al.^^ obtained at 50 K and downloaded 
as GIF file from the IGSD,^^ with a* = 6* = 5.507 A 
and c* = 7.796 A. The starting rotation angle for TiOe 
octahedra was 2.1° while the c/a — 1 = 10 x 10~^. 
Please note the geometries were only starting points; as 
mentioned above all geometries were optimized with the 
method/basis set under consideration. In order to avoid 
introducing any errors coming from size effects or k-space 
integration, the calculated properties of the AFD super- 
cell are always compared with a 20 atoms supercell con- 
structed from four cubic primitive cells (without octa- 
hedral rotation or tetragonality) fully relaxed using the 



same k -point mesh. It should be noted that the supercell 
in the cubic phase is a local minimum and is higher in en- 
ergy than the supercell in the AFD phase for all reported 
calculations. The final (reported) values were deter- 
mined from Ti-02-Ti angle measurements, and any octa- 
hedral tilts can be estimated by measuring the Ti-Oi-Ti 
angles (O^'s subscript was defined in Figure 1). Finally, 
all geometric visualization was done using Gauss View. 

III. BASIS SET EFFICIENCY FOR SrTiOs 

The challenge in selecting a basis set is always bal- 
ancing accuracy with computational cost. In molecular 
calculations, the computational cost of a gaussian basis 
set is determined by the number of functions used, while 
in PBG calculations the spatial extent or diffuseness of 
the basis set also plays a major role. The more diffuse a 
basis set is, the larger the chunk of matter that must be 
included in the calculations to avoid numerical issues. 

Goupled with the argument that the long density tail is 
more necessary for molecular work than work in extended 
systems, it becomes obvious that basis sets developed 
for non-periodic calculations can require modification for 
PBG use. This section describes the basis set optimiza- 
tion/modification procedure we employed to find the ap- 
propriate Gaussian basis sets to simulate periodic STO 
while keeping within reasonable computational expense. 
We based our evaluations of a basis set's accuracy on cu- 
bic STO results using the Heyd-Scuseria-Ernzerhof^'''^^ 
screened Goulomb hybrid density functional (HSEOG).^'' 

The obvious starting point was the basis sets used in 
previous calculations/studies of bulk STO, including: 

• Gaussian-type basis sets published by Piskunov 
et al.^^ in 2000, optimized using the Hartree-Fock 
(HF) and density functional theory (DFT) with 
Hay-Wadt pseudopotentials^^"^^ for Sr and Ti, 
denoted here as PI. 

• The subsequently improved version of PI pub- 
lished by Piskunov et al.^^ in 2004, which expands 
of PI by adding polarization t/- functions to oxygen 
and making the Ti s and p functions more diffuse, 
denoted here as P2. 

Tests on PI and P2 were done with HSE06, because 
it has been found to give the best results versus exper- 
iment for both structural and electronic properties'^ in 
older calculations. Both PI and P2 reproduce the ex- 
perimental equilibrium lattice constants^^ (see Table I) 
almost perfectly. Gubic STO modeled with PI has a 
slightly higher bulk modulus compared to P2, although 
the difference between the two basis sets in fairly mini- 
mal for structural properties. A more important effect is 
observed for the electronic properties: PI and P2 over- 
estimate the direct band gap of STO by 0.12 and 0.05 
eV respectively, and seriously overestimate the indirect 
band gap by 0.28 and 0.21 eV. 
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TABLE 1. The electronic and structural properties of cubic 
SrTiOs computed with HSE06^^ and different basis sets. 
Please see the text for basis set naming conventions. 

Basis set PI P2 SZVP TZVP Experiment 



Direct gap (eV) 3.87 3.80 3.59 3.59 3.75^ 
Indirect gap (eV) 3.53 3.46 3.18 3.20 3.25^ 
ao(A) 3.900 3.908 3.887 3.902 3.890^, 3.900^ 

B(GPa) 198 194 204 193 179^, 

179zb4.6^ 



^ Reference 63. 
^ Reference 56. 
^ Reference 64. 
^ Reference 65. 

It is easy to see that the P2 basis set employed with 
HSE06 lead to results that are closer to experiment than 
PI, a fact noted by Piskunov^^ for a number of func- 
tionals. The more important point is that increasing the 
size/quality of the basis set made a noticeable change in 
the results; the immediate question is whether another 
increase in basis set size would bring about similar im- 
provement. In other words, using polarization c?-orbitals 
for O and diffuse functions for Ti improved the HSE06 
results, and imply that further improvement could poten- 
tially be achieved if more basis set improvements are im- 
plements, e.g. including titanium core electrons and/or 
adding more diffuse functions for oxygen. 

We decided to optimize some of the Def2-^^ series of 
Gaussian basis sets for use in bulk STO calculations. The 
original Def2- basis sets for the atoms of interest in this 
project included small- exponent diffuse functions {amin 
less than 0.10) that are spatially quite extended; as stated 
above, this long tail is necessary to improve the DFT re- 
sults for molecules but not necessary for crystals.^ '''^^ Ba- 
sis sets with large spatial ranges dramatically slow down 
the calculation of Coulomb contributions to the total en- 
ergy of crystals. Thus, to be useful in PBC calculations, 
Def2- basis sets must be modified by removing some of 
the most diffuse functions. 

The series of Def2- basis sets are available up 
to quadruple zeta valence quality for a large set of 
elements. ^^'^^ In the original optimizations, the oxygen, 
strontium and titanium basis sets were optimized (us- 
ing HE and DET) versus the properties of SrO, TiO and 
Ti02 molecules. Strontium has the inner shell electrons 
replaced with small core pseudopotentials,^^ while the 
other two atoms utilize all electron basis sets; this differs 
from PI and P2 which uses pseudopotentials on tita- 
nium as well. In general, Def2- basis sets are larger and 
more expensive than PI and P2 basis sets, but are ex- 
pected to give a better representation of both phases of 
STO due to greater "completeness." 

To make a Def2- basis set applicable to PBC, the 
first step is selecting a maximum allowable diffuseness, 
or equivalently the smallest acceptable Gaussian orbital 



exponent, Ofmin- The larger the value of a^nin^ the faster 
the calculations become, but if amin is set too high, sig- 
nificant degradation of physical property prediction re- 
sults. After the threshold is defined, one pass is made 
through the basis set to reset all a < amin to a^in^ and 
then a second pass is made through the basis set to re- 
move any redundancies. Note that after modifying or 
deleting an element of a contracted basis set, we rely on 
the code's internal renormalization code, i.e. no attempt 
is made to reoptimize contraction coefficients. 

We first began with the largest Def2- basis sets, Def2- 
QZVP and Def2-QZVPP, but these were found to be 
computationally intractable for bulk STO even for amin 
as big as 0.2, and previous experience have shown that 
amin larger than 0.2 causes physically unacceptable re- 
sults. We then moved to the smaller basis sets, Def2- 
TZVP and Def2-SZVP. We first set amin= 0.12, but 
found this made the calculations very slow. Our tests 
showed that amin= 0-15 constitutes a more computation- 
ally efficient choice without important loss in accuracy. 

Henceforth, the Def2-TZVP and Def2-SZVP, with 
amin modified and redundant s functions removed, will 
be denoted TZVP and SZVP, respectively. 

Table I summarizes the calculated electronic and struc- 
tural properties of cubic STO using our basis set modi- 
fications as well as the aforementioned PI and P2. The 
optimized basis sets SZVP and TZVP give an over- 
all excellent agreement with experiment direct band 
gaps are now underestimated by 0.16 eV while indirect 
band gaps are now underestimated by ^0.05 eV. These 
two new basis sets are larger than the previously utilized 
PI and P2, are more accurate for indirect gaps, as well 
for other measured properties, and due to their greater 
size are expected to be closer to the upper limit of HSE06 
accuracy for this system. Note also that the electronic 
properties of STO remain almost unchanged by moving 
from a SZVP to TZVP basis set. The deviation from 
the experimental lattice constant do not exceed 0.07 and 
0.3% for SZVP and TZVP respectively but is more sub- 
stantial for the bulk modulus reaching 14% for SZVP 
and 8% for TZVP. Einally, the same series of basis set 
optimizations were also performed using HISS and M06L 
functionals, which lead to the same conclusions regarding 
the basis set efficiency; these are not presented here for 
space reasons. 

Before moving on to the results section, a brief men- 
tion of the expense of the various basis sets should be 
included. In term of relative CPU time, one SCE cycle 
takes about 12 units for TZVP compared to 6 units for 
SZVP and 1 unit for P2. All of these basis sets still 
have potential uses; SZVP or P2, for example, might 
be very useful for a rapid investigation of the electronic 
properties of some complex STO systems. But, in term of 
completeness, TZVP is the most complete and the clos- 
est to the planewave basis set limit, followed by SZVP, 
then P2. 
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IV. RESULTS: BASIS SET AND FUNCTIONAL 
EVALUATION 

In this section we present the calculated properties of 
SrTiOs, always discussing the results of each functional 
using the TZVP basis set first, followed with a discussion 
of the sensitivity of the functionals to smaller basis sets, 
namely SZVP and P2. 



A. Structural properties of cubic SrTiOs 

The calculated equilibrium lattice constants and the 
bulk moduli of cubic STO using different functionals and 
basis sets are reported in Table 11. Unless specified, the 
deviation of theory from experiment will be always re- 
ferred to the data of Abramov et al,^^ i.e. treated as 
the target value. Focusing first on the TZVP results, 
we observe that the screened hybrids HSE06 and HISS 
give lattice parameters in excellent agreement with ex- 
periment. 

The calculated bulk modulus using HSE06 is fairly 
close to the experimentally reported values, although 
overestimated by 8%. (The same magnitude of overes- 
timation have been also reported in the HSE planewave 
calculations of Wahl et al.^^.) However, a larger bulk 
modulus overestimation of 15% is observed for HISS, 
which constitutes the largest deviation from experiment 
among all the studied functionals. 

M06L and revTPSS predict slightly higher equilibrium 
lattice constants than screened hybrids do, but their bulk 
moduli are closer to experiment, with revTPSS being 
especially close. LSDA underestimates the lattice con- 
stant by 0.03 A, while PBE predicts lattice constants 
0.05 A larger than experiment. PBEsol is in excellent 
agreement with the experimental lattice constant. Thus 
PBEsol corrects the LSDA underestimation and the PBE 
overcorrection to LSDA for lattice constants; in addition, 
the PBEsol bulk modulus deviate by less than 3% from 
experiment, while LSDA and PBE are off by 11% and 
12%, respectively. This is an example of PBEsol meet- 
ing its purpose, as it improves the PBE lattice constant 
and bulk modulus for the cubic phase, approaching very 
closely the experimental data. 

Turning now to the functional sensitivity to basis set 
size, we observe from the HSE06 results that SZVP ba- 
sis set predict bond lengths that are very slightly shorter 
than the TZVP and a bulk modulus that is 6% higher. 
As such, SZVP predicts SrTiOa to be 14% harder than 
experiment. P2 behave in the opposite direction, pre- 
dicting slightly longer bonds when compared to TZVP, 
while the bulk moduli are only IGPa higher. From ta- 
ble II, this sensitivity of HSE06 to smaller basis set can 
be generalized to M06L, revTPSS and the semilocal func- 
tionals LSDA, PBE and PBEsol. 

Finally, it should be noted that PBEsol results offer 
the best agreement with experimental structural proper- 
ties^^ of SrTiOs among all the studied functionals with 



the TZVP basis set, followed by the screened hybrid 
HSE06 and the meta-GGA revTPSS. 



B. Electronic properties of cubic SrTiOs 

The computed electronic properties of SrTiOs are sum- 
marized in table HI. As expected, HSE06 gives an excel- 
lent estimate of the electronic properties when used with 
the large TZVP basis sets. Deviations from the experi- 
mental values are 0.15 eV for the direct gap and 0.05 for 
the indirect gap. A cursory glance over the rest of ta- 
ble HI indicates that no other functional was comparable 
to HSEOG's efficacy for band gaps, i.e. everything else we 
tried had much larger errors. 

The midle range screened hybrid HISS tend to over- 
estimate the direct and indirect band gaps by 0.35 and 
0.73 eV, respectively. M06L and revTPSS tend to un- 
derestimate both band gaps, by an average of ~1.2 and 
~1.4 eV respectively. The semilocal functionals LSDA, 
PBEsol, PBE underestimate the experimental band gaps 
by an average of 45% or 1.5 eV. This was expected, and 
is in agreement with the behavior observed earlier in the 
literature for this system. It can be easily seen from 
these results that HSE06 is the best functional choice for 
investigating this system. 

Turning to basis set sensitivity, it can be seen from the 
HSE06 numbers that band gaps are nearly unaffected by 
using the smaller SZVP basis sets, but when used with 
the still smaller P2 basis set, direct and indirect band 
gaps increase by ^0.25 eV versus TZVP. The predicted 
direct band gap becomes closer to experiment when using 
P2 and HSE06, probably due to a cancellation of errors 
effect, while the indirect band gap is noticeably worse. 
This same sensitivity holds for almost every other func- 
tional, with SZVP and TZVP giving about the same 
band gaps and P2 opening the band gaps up by a few 
tenths of an eV. M06L appears to be slightly less sensi- 
tive; no obvious reason for this exists. 



C. Stability of the AFD phase of STO 

This section examines the stability of the AFD phase 
of STO as calculated by the various functionals and ba- 
sis sets previously tested for the cubic phase. The func- 
tional/basis set combinations tested face the challenge of 
predicting the AFD octahedral rotation angle, 0, as well 
as the tetragonality parameter c/a, which as shown in 
section I is not trivial. The AFD phase order parame- 
ters are evaluated from the relaxed 20-atoms AFD super- 
cells as described in section II. The performance of each 
functional with TZVP, and then an analysis of the func- 
tional's sensitivity to the smaller basis sets is presented 
in turn. 

Figure 2 shows that the screened hybrid functional 
HSE06 is excellent for the structural properties of AFD, 
as it was for the cubic phase. Both the rotation angle 
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TABLE 11. Computed lattice parameter ao(A) and bulk modulus B(GPa) for Cubic STO using different combinations of 
functionals and basis sets compared to experiment. 

HSE06 HISS M06L revTPSS LSDA PEE PBEsol Experiment 



ao(A) 3.890^, 3.900^ 

TZVP 3.902 3.883 3.925 3.921 3.862 3.941 3.897 

SZVP 3.887 3.869 3.909 3.903 3.845 3.924 3.881 

P2 3.908 3.891 3.930 3.920 3.870 3.946 3.903 

B(GPa) 179^, 

TZVP 193 206 187 180 201 169 184 179=b4.6^ 

SZVP 204 218 198 193 214 180 196 

P2 194 205 191 184 203 173 187 



^ Reference 56. 
^ Reference 64. 
^ Reference 65. 



TABLE IIL Direct and indirect band gaps computed for Cubic SrTiOs using different basis sets and functionals compared to 
experiment. 



HSE06 


HISS 


M06L 


revTPSS 


LSDA 


PBE 


PBEsol 


Experiment 



Direct gap(eV) 3.75" 

TZVP 3.59 4.39 2.51 2.24 2.08 2.11 2.10 

SZVP 3.59 4.45 2.53 2.28 2.12 2.14 2.14 

P2 3.80 4.56 2.63 2.52 2.34 2.33 2.34 

Indirect gap(eV) 3.25" 

TZVP 3.20 3.98 2.09 1.87 1.75 1.74 1.75 

SZVP 3.18 4.03 2.10 1.89 1.76 1.75 1.76 

P2 3.46 4.22 2.24 2.17 2.04 1.99 2.02 



and the c/a ratio are in very good agreement with ex- 
periment. These properties are not significantly affected 
when SZVP is used, but HSE06/P2 predicts a very very 
small angle for the AFD phase, while retaining a good 
c/a. This is one area where TZVP noticeably outper- 
forms P2 with HSE06. 

HISS and revTPSS behave as HSE06 for both TZVP 
and SZVP, giving a good estimate of both order param- 
eters. However, they demonstrate a higher sensitivity 
to the smaller P2 basis set and required the use of a 
very stringent convergence criterium to finally relax the 
structure back to a pseudocubic phase with ?=^0. On 
the other hand, M06L predicts the AFD phase to be un- 
stable, and relaxes to a non-rotated structure regardless 
of the basis set used. 

The semilocal functionals LSDA, PBEsol and PBE aU 
overestimate the tetragonality of the AFD phase by pre- 
dicting and c/a almost twice the size of the experimen- 
tal results. The highest overestimation was observed for 
LSDA, followed by PBEsol then PBE. Note that our re- 
sult here is in excellent qualitative agreement with the 
behavior found in the planewave calculations of Wahl 
et al;^^ quantitatively, however, the LSDA, PBEsol and 



PBE octahedral angles with TZVP are 25-30% lower 
than the planewave results^^^^^ (for a detailed numer- 
ical comparison see table IV C, and ref. 23 has addi- 
tional comparison with experiment). Similar behavior 
has been recently published''^ for LSDA calculation with 
finite-range numerical atomic orbitals using a double-^ 
polarized basis set. This indicates that localized basis 
sets tend to reduce the AFD octahedral rotation com- 
pared to plane-waves but do not succeed to suppress the 
DFT overestimation. 

When used with the SZVP basis sets, the LSDA, PBE 
and PBEsol rotation angles are larger than the TZVP 
ones. Furthermore, when LSDA, PBEsol and PBE are 
used with the P2 basis set, we observe a small and co- 
herent reduction in the octahedral rotation angle of the 
AFD structure compared to TZVP results. This demon- 
strates that semilocal functionals have different degrees 
of sensitivity to the quality of the localized basis sets 
used, but the error caused by functional choice is always 
the more important source of error. Thus the functionals 
examined here will lead to exaggerated AFD values for 
all basis sets considered. 
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^■TZVP 
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^■P2 
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FIG. 2. (Color online) Performance of different functional/basis set combinations in predicting the order parameters of the 
AFD phase transition in STO. Dashed lines depict the experimental octahedral angle measured at 4.2 K from Ref. 11 (left) 
and the tetragonality parameter obtained at 50 K from Ref. 12 (right). 



DISCUSSION: PHYSICAL PROPERTIES OF 
STO 



Before talking about specific issues, there are a few 
general conclusions we can reach from examining the re- 
sults in section IV 



1. HSE06/P2 did a good job in describing accurately 
the structural properties for the cubic phase as 
well as providing a descent estimation of the band 
gap. However, the failure of HSE06/P2 to correctly 
model the structure of the AFD phase indicates 
that it must be abandoned as a useful combination 
for this and related systems. 



2. HSE06/SZVP has the drawback of predicting a 
stiffer SrTiOa in the cubic phase, although it pre- 
dicts electronic properties as well as TVZP. It also 
predicts a stiffer AFD structure, but the octahedral 
angle and cja parameters are very good. 



3. HSE06/TZVP gave the best agreement with ex- 
periment for the cubic phase and for the AFD 
phase. It is definitely the most reliable combina- 
tion of functional and basis set among all studied 
variations. Thus HSE06/TZVP can be used with 
confidence on more complicated structures, as well 
as to understand the change in the electronic struc- 
ture during the cubic to AFD transition for this 
system. More concisely, we believe that this com- 
bination is an accurate enough functional in a good 
enough basis set to explain phenomena in metal ox- 
ides. 




FIG. 3. Band structure of cubic SrTiOs calculated with 
HSE06/TZVP. The dashed line depict the Fermi level lying 
at the valence band maximum (R special point.) 




FIG. 4. (Color online) Total electronic density of states 
(DOS) of cubic SrTiOs calculated with HSE06/TZVP. Pro- 
jected density of states (PDOS) of the main contributing or- 
bitals are also shown. 
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TABLE IV. Structural and electronic properties of the antiferrodistorsive phase of SrTiOs compared to previously simulated 
data and experiments, a* and c* are the lattice parameters of the AFD supercell and c/a = c* /y/2a*. AE = Ecubic — Eafd 
represent the gain in total energy after the cubic to AFD phase transition, while AEg denote the corresponding increase in the 
band gap. 





LSDA 


PBE 


PBEsol 


HSE06 


HISS 


revTPSS 


M06L 


Experiment 


a* (A) 
Present 
Ref. 23^ 


5.449 
5.440 


5.568 
5.562 


5.500 
5.495 


5.515 
5.515 


5.448 


5.543 


5.551 


5.507^ 


c*(A) 

Present 
Ref 23^ 


7.727 
7.755 


7.900 
7.897 


7.812 
7.818 


7.809 
7.808 


7.772 


7.846 


7.862 


7.796^ 


(c/a - 1) X 10""^ 
Present 
Ref. 23^ 
Others 


27 
80 
40^ 


32 
40 


44 
60 


12 
10 


14 


7.6 


7 


10^ 


en 

Present 
Ref. 23^ 

Ot.bprft 


4.14 
6.05 

8 40^ 6"^ 


3.54 
4.74 


3.81 
5.31 


2.01 
2.63 


1.92 


2.01 





2.01zb0.07^ 
2.1^ 


AE X 10"^ (eV) 
Ref. 23^ 


1796 
1900 


854 
700 


44 
1100 


35 
200 


578 


258 


122 




Indirect band gap (eV) 

Present 

Ref. 23^ 


1.820 
1.970 


1.787 
1.790 


1.808 
1.930 


3.227 
3.110 


3.995 


1.890 


2.060 


3.246^^ 
3.160^^ 


AEg (meV) 
Present 
Ref. 23^ 


75 
160 


49 
10 


58 
110 


27 
40 


15 


15 


30 


50^ 



^ plane- wave calculation using a different HSE screening parameter. 
^ Reference 25. 
^ Reference 24. 

^ Reference 73 using numerical atomic orbit als. 
^ Reference 12 (at 50K). 
^ Reference 11 (at 4.2K). 

s Reference 71 difference between 85 K and 8 K measured gaps. 



A. Band structure alteration by the AFD phase 
transition 

The band structure of the cubic unit cell of STO com- 
puted with HSE06/TZVP is shown in figure 3, with 
the high symmetry points F = (0,0,0), X = (0,^,0), 
M = (^, ^,0) and = (^, ^, ^) labeled, in the first Bril- 
louin zone of the simple cubic system. The dashed line 
depict the Fermi level lying at the valence band maxi- 
mum at the R point. 

Our band structure agrees qualitatively with previous 
band structures from LSDA/PW calculations, which can 
be seen (for example) in Fig. 5 of Ref. 25, as well as the 
B3PW/P2 band structure in Ref. 33 Fig. 2(a), with 
the exception of a few details. Our direct band gap 



(F F) of 3.59 eV and indirect gap {R F) of 3.2 eV 
are in better agreement with experiments^ compared to 
the underestimation observed in the LSDA/PW gaps 
and the overestimation found with B3PW/P2. Thus 
for a DFT approach, this diagram is the best band 
structure to date. (An even more accurate band struc- 
ture was computed using the experimental lattice con- 
stant of SrTiOs by mean of post-LSDA quasiparticle self- 
consistent GW (QSGW) correction to the band structure 
by Hamann and Vanderbilt.''^) 

Figure 4 shows the total density of states (DOS) as 
well as the projected density of states (PDOS) on every 
atomic orbital. The PDOS of oxygen represents the sum 
of the contributions of all three oxygen atoms in the cu- 
bic unit cell. In the energy window shown here, the DOS 
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is dominated by oxygen 2p, titanium 3d and strontium 
Ad states. (All the remaining orbitals have a negligible 
contribution, so their PDOS are not shown.) The valence 
band (VB) from to —6 eV is dominated by oxygen 2p 
states, with a small contribution from titanium 3d states 
in the range —3 to —6 eV. The conduction band (CB) 
is clearly dominated by titanium 3d in the energy range 
3.2—7 eV, with a smaller contribution coming from the 
3 oxygen 2p states as well. The admixture in the VB 
and CB between the titanium 3d and oxygen 2p orbitals 
demonstrates that the Ti— O bonds have a partially cova- 
lent character with a small degree of hybridization. (This 
behavior has been noted in previously published data.^^) 
Between 7—9 eV, the spectrum is the sum of contribu- 
tions from oxygen 2p, titanium 3d and strontium M or- 
bitals. The higher energy region in the CB (9—12 eV) is 
dominated by strontium Ad orbitals with small contribu- 
tions from titanium 3d and with oxygen 2p vanishing at 
around 10.5 eV. 

Figure 5 compares the total electronic densities of 
states for the cubic and AFD supercells. As a general 
trend, the cubic to AFD phase transition does not lead to 
a significant modification in the total DOS; both the va- 
lence and the conduction bands experience a slight shift 
to higher energies (horizontal arrows) together with some 
small modifications, indicated by vertical black arrows 
in Figure 5. However, the VB shift does not affect the 
peak at the VB maximum, while a very small shift to 
higher energies is observed for the CB minimum, indi- 
cating that the band gap increase by ^27 meV after the 
transition. This same behavior holds for all the function- 
als/basis sets combinations tested, and is in line with 
some experimental observations^^ reporting very small 
changes in their measured band gaps due to the cubic- 
to-tetragonal structural transition. Further confirmation 
of this physical effect can be seen in recent photolumines- 
cence measurements,^^ which reported that the band gap 
increased by 50 meV when temperature decreased from 
85 K to 8 K, which is a temperature range over which the 
AFD rotation would go from incomplete to nearly total. 

A more detailed comparison between the PDOS for 
each atomic orbital can give a better understanding of 
the origin of these modifications. It is important to men- 
tion that in the AFD supercell, there is one nonrotating 
Oi atom and two rotating O2 oxygens for every Sr and 
Ti atom. Concentrating on the oxygen 2p orbitals, we 
observe that the non-rotating Oi atoms are nearly un- 
changed in the PDOS compared to the cubic phase, with 
the exception of a tiny shift to higher energy (see the inset 
in Figure 6), which can be attributed to the elongation of 
the cell along the z-axis. However, the O2 demonstrate a 
much more significant shift to higher energies, along with 
changes in the height and width of some peaks. This is 
mainly caused by the octahedral rotation involving O2 
atoms. The titanium 3d and strontium M spectra expe- 
rience the same aforementioned shift to higher energies 
in the VB and the CB due to the elongation of the lat- 
tice, with a few noticeable changes in the titanium 3d 
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FIG. 5. (Color online) Modification in the total DOS of STO 
upon the cubic to AFD phase transition. Horizontal arrows 
indicate the direction of the energy shift, while vertical arrows 
point to the most important changes. 



spectrum at —2.9 as well as between 5 and 6.5 eV. Most 
of the modifications observed in the total DOS, with the 
exception of few originate from the changes in the O2 2p 
and Ti 3d spectra with the O2 being far more important. 



B. The effect of the HSE screening parameter, uj 

Relying on the assumption that plane- waves are much 
closer to the infinite basis set limit than the Gaus- 
sian basis sets we used, it is useful to compare our 
HSE06/TZVP results with the HSE plane-wave results. 
To our knowledge, only Wahl et al.^^ have published data 
using plane-waves and a Heyd-Scuseria-Ernzerhof^^'^^'^^ 
style screened Coulomb hybrid density functional for this 
system. However, a direct comparison with our present 
data is not possible because Wahl et al. used a different 
screening parameter in their calculations. 

Briefiy, the HSE functional partitions the coulomb po- 
tential into short-range (SR) and long-range (LR) com- 
ponents: 

(1) 
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FIG. 6. (Color online) Modification of the partial densities of states (PDOS) for O 2p, Ti 3d and Sr 4d Left: Valence band 
Right: Conduction band 



The screening parameter uj defines the separation 
range, as it controls the distance at which the long-range 
nonlocal interaction becomes negligible, i.e. it "turns 
off" exact exchange after a specified distance. Wahl 
et al. used cji =0.159 a.u.~^ ^ effectively using an HSE- 
style functional, but not either of the functionals HSE03 
or HSE06. Krukau et al.^^ applied HSE while varying 
0.11< UJ < 0.20 to a number of bulk metals and semi- 
conductors. They concluded that a small increase of uj 
substantially lowers the calculated band gaps and the 
smaller value omega takes in this range, the closer the 
calculated band gaps and lattice constants are to the ex- 
periment. Based on the above, Krukau et al.^^^ recom- 
mended co'2=0.11 a.u.~^ for both the HF and PBE part 
of the exchange. This is the value we used in all our cal- 
culations, and this value is part of the formal definition 
of HSE06. However, in order to make a comparison be- 
tween our HSE(cj2)/TZVP and the HSE(a;i)/PW data of 
Wahl et al., we must perform a HSE/TZVP calculation 
with uji and isolate the screening parameter effect on the 
calculated properties of SrTiOa. 

Table V shows that the HSE(cJi)/TZVP lattice con- 
stant and bulk modulus changes very slightly by decreas- 
ing the screening parameter from uji to UJ2'- the change 
is 0.001 A and 1 GPa respectively. A much more sig- 
nificant effect is, however, observed for the band gaps: 
decreasing the screening parameter by 50% { uoi ^ UJ2) ^ 



TABLE V. Variation of the cubic STO lattice parameter(ao 
in A ), bulk modulus (B in GPa), and direct {Eg) and indi- 
rect {Eg) band gaps (in eV) by decreasing the HSE screening 
parameter from cji=0.159 a.u.~^ to CJ2=0-11 a.u.~^ . Results 
are from our Gaussian basis set (TZVP) and the plane-wave 
(PW) calculations in Ref. 23. 





Gaussian 




PW 


Experiment 


LOi — 


(J02 




— > (J02 


ao 


3.903 


3.902 


3.904 


3.903^ 


3.890^, 3.900^ 


B 


192 


193 


192 


193" 


179^, 179=b4.6^ 




3.37 


3.59 


3.47 


3.67" 


3.75^ 


K 


2.96 


3.20 


3.07 


3.27" 


3.25^ 



^ Reference 56. 
^ Reference 64. 
^ Reference 65. 
^ Reference 63. 

^ Estimated values if 6^2=0.11 a.u~^ was used in plane- wave 
calculations of Ref. 23. 



lead to an increase in the band gaps, effectively a rigid 
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shift of 0.22 and 0.24 eV for the direct and indirect band 
gaps respectively. If examined from the other direction, 
decreasing the screening parameter from uoi to UO2 (with 
HSE/TZVP) tends to bring the band gaps closer to the 
experiment (see table V), which suggests that UJ2 pro- 
vides better agreement with experiment than uji does. 
The same structural changes and band gap shifts were 
also found for the small basis sets SZVP and P2, which 
are not presented here and which demonstrate that this 
effect is completely independent from the basis set used. 
Finally, the HSE(cj2)/TZVP band gaps are very close to 
the HSE(cj2)/PW values we estimated, suggesting that 
our TZVP basis set is very close in quality to the pre- 
viously used plane waves, and thus is closer to the basis 
set limit. 

This section is contains one of the most important re- 
sults of this paper, and as such should be clearly restated. 
If we use the same version of HSE used in plane wave 
studies, we can show that our TZVP is a high quality 
basis set as it matches the excellent basis set plane wave 
results. If we use the proper uj in HSE with our basis 
set, we arrive at the best results/smallest errors versus 
experiment ever reported for SrTiOs. 

Finally, it should be noted that this is not an ad hoc 
parameterization of (jO to give the best results for this 
study. We were able to obtain results that closely match 
experiment by using a demonstrably high quality basis 
set and a parameter in the density functional determined 
by a large test bed of structures and properties. 



C. Screened hybrids compared to regular hybrids. 



and B3LYP, done with the P2 basis set. There are no- 
ticeable differences between the results of HSE06 and the 
regular hybrids, with HSE06/TZVP giving band gaps 
very close to experiment while regular hybrids used with 
P2 overestimate the gap, especially the indirect band 
gap. The band gap overestimation is of same magni- 
tude we observed in section IV B for HSE06/P2 as well 
as all the other functionals tested on STO with P2. 
This suggests that P2 is also behind the band gap over- 
estimation in the regular hybrids data reported in the 
literature.^O'^^'^^ By comparing the P2 and TZVP band 
gaps from table III, we can deduce that the P2 basis set 
has an effect (versus a large basis set) of increasing the 
direct and indirect band gaps by average values of 0.22 
and 0.28 eV respectively. By applying this P2^TZVP 
basis set correction to the regular hybrid B3PW/P2 and 
B3LYP/P2 band gaps (see the corrected values in ta- 
ble VI), band gaps are brought closer to the experimen- 
tal values, and thus closer to the HSE06/TZVP results 
as well. Consequently, differences in the computed elec- 
tronic properties of HSE06 and B3PW and B3LYP are 
considerably attenuated and suggest that the screened 
hybrid HSE06 is comparable in accuracy with regular hy- 
brids for STO, while being much more computationally 
efficient. 

The final issue to discuss is the comparison of the 
structural and elastic properties of STO computed with 
HSE06 versus regular hybrids. Perovskite crystals in the 
cubic structure have only three independent elastic con- 
stants, namely Cn, C12 and C44, as well as a bulk mod- 
ulus: 



5 = i(Cn+2Ci2) 



(2) 



TABLE VL Our most converged direct {Eg) and indirect {El) 
band gaps (in eV) for cubic STO alongside previously pub- 
lished hybrid functional results done with the P2 basis set. 
Regular hybrids data are corrected according to the basis set 
sensitivity effect deduced in section IV B. 



functional /basis 



E": 



Ori. Corr.^^'^"^^'^^ Ori. Corr.^^'^"^^'^^ 



Exp.^3 


3.75 


3.25 


Present HSE06/TZVP 


3.59 


3.20 


Ref. 33 B3PW/P2 


3.96 3.74 


3.63 3.35 


B3LYP/P2 


3.89 3.67 


3.57 3.30 


Ref. 75 B3PW/P2 


4.02 3.80 


3.70 3.42 


Ref. 20 B3PW/P2 




3.63 3.35 


Ref. 76 B1-WC/P2^ 


3.91 


3.57 



^ P2 basis set with all electrons for Ti, basis set correction 
cannot be applied. 



Table VI summarizes the calculated band gaps of 
HSE06/TZVP and compares them with previously pub- 
lished gaps computed with the regular hybrids B3PW 



We calculated the elastic constants of STO using 
HSE06/TVZP, following the methodology described in 
Ref. 77. Ideally we would like to compare our cubic 
elastic constants calculated at K with low temperature 
data, but as experimentally the cubic structure turns to 
a tetragonal structure below a transition temperature, 
making any comparison of this kind impossible. Exper- 
imentally, Bell and Rupprecht''^ found that the elastic 
constants of STO measured between 303 and 112 K obey 
the following empirical relations: 





= 334. 1[1 - 


2.62 X 10" 




-Ta) 


C12 


= 104.9[1 - 


1.23 X 10" 


-\T 


-Ta) 


C44 


= 126.7[1 - 


1.30 X 10" 


-4(T 


-Ta) 



0.0992 

{T - To) 
0.1064 

(T - Ta) 

0.1242 

{T - Ta) 



where the elastic constants are in GPa, T is the temper- 
ature and Ta=108 K is the critical temperature. On and 
C44 reach their maximum values at 133 K where STO is 
still cubic, then they start to decrease as — 1/(T — T^) 
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in the region around the transition temperature, in con- 
trast, Ci2 continue to increase as 1/(T — T^) in the same 
temperature range. 

Since we don't know at which temperature the change 
from cubic to tetragonal phase begins to take place, it 
is better to limit our comparison with data measured 
at 133 K and above. Table VII summarizes our results 
and compares them with experiment as well as previ- 
ously published results with B3PW/P2 and B3LYP/P2. 
HSE06/TZVP provides excellent lattice constants but 
predicts the bulk modulus to be 8% higher than exper- 
iment. The elastic constants from HSE06/TZVP over- 
estimate the experimental data at room temperature by 
10% and the 133 K data by 6%; this was expected given 
the overestimation of the bulk modulus. The B3PW hy- 
brid also gave very good lattice constant and bulk mod- 
ulus, but the calculated elastic constants are lower than 
the room and the low temperature experimental values. 
B3LYP predicted a lattice constant higher by 1%, a good 
bulk modulus, and offers the best agreement with the 
low temperatures elastic constants. In summary, none 
of the screened or regular hybrids considered was able 
to give simultaneously excellent bulk moduli and elastic 
constants, still HSE06/TZVP offer the best compromise 
between efficiency, accuracy and speed. 

TABLE VIL Calculated elastic constants with HSE06/TZVP 
for cubic STO compared to experiment and previously pub- 
lished results with the regular hybrid functional B3PW and 
the P2 basis set from Ref. 33. ao is in A, B, Cn, C12 and 
C44 are in GPa 





ao 


B 


Cii 


C12 


C44 


HSE06/TZVP 


3.902 


193 


351.4 


113 


137.3 


B3PW/P2 


3.900 


177 


316 


92.7 


120.1 


B3LYP/P2 


3.940 


177 


328.3 


105.7 


124.6 


Exp. 


3.890^^ 


179^^ 


317.2 


102.5 


123.5^ 




3.900^^ 


179zb4.6^^ 


330 


105 


126^ 




3.910 


184 






128"^ 



^ Ref. 78 at room temperature. 

^ Ref. 78: max. measured values for Cn and C44 at 133 K, C12 

increase further as temperature drop. 
^ LandoltBornstein Group III Condensed Matter 2002 vol 36, 

subvol V (Berlin: Springer) chapter lA (Simple 

Perovskyte-Type Oxides) pp 11647 



VI. CONCLUSION 

We used the ab-initio code GAUSSIAN to simulate the 
properties of SrTiOa (STO) using a large spectrum of 



functionals, from LSDA, GGAs (PBE and PBEsol) and 
meta-GGAs (M06L and revTPSS) to modern range- 
separated hybrid functionals (HSE06 and HISS); assess- 
ing their ability in predicting the properties of the cubic 
and the AFD phases of STO. 

We found that pure DFT functionals tend to overesti- 
mate the octahedral rotation angles of the AFD phase, in 
agreement with previously reported results in the litera- 
ture using plane- wave basis sets of comparable quality. 
Also, basis sets of low quality tend to inhibit the tetrag- 
onality of the AFD phase and sometimes even suppress 
it, regardless of the functional used. We therefore con- 
structed a localized basis set of sufficient completeness (or 
size) to correctly simulate the TiOe octahedral rotation 
and the cubic phases of STO. We also evaluated the band 
gap errors arising from the use P2 basis set and from the 
magnitude of the HSE screening parameter oj. By apply- 
ing our basis set and uu corrections to the previously pub- 
lished work with regular and screened hybrid functionals 
on STO, we showed that the discrepancies between pub- 
lished simulated data can be explained and that hybrid 
functionals used with sufficiently big Gaussian-type basis 
sets can give results comparable with plane-wave calcu- 
lations and in excellent agreement with experiment. 

The screened hybrid functional HSE06 predicts the 
electronic and structural properties of the cubic and 
AFD phase in very good agreement with experiment, 
especially if used with high quality basis set TZVP. 
HSE06/TZVP is the most reliable combination of func- 
tional and Gaussian basis set for STO which is computa- 
tionally tractable with the current computer power. It is 
accurate enough to enable us to understand the changes 
in the band structure during the cubic to AFD phase 
transition. The success of HSE06/TZVP encourages its 
use on more complicated cases like the bond breaking 
and over binding and defect formation, where the basis 
set completeness is expected to play a major role. 
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